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ABSTRACT 



We perform calorimetry on the bright gamma ray burst (GRB) GRB 030329 
by fitting simultaneously the broadband radio afterglow and the observed after- 
glow image size to a semi-analytic magnetohydrodynamical (MHD) and afterglow 
emission model. Our semi-analytic method is valid at all relativistic regimes and 
is successfully used to fit archival measurements of the afterglow flux from 1 
day to 5.5 years after the burst. Values for the initial burst parameters are de- 
termined and the nature of the circumburst medium is explored. Additionally, 
direct measurements of the lateral expansion rate of the radio afterglow image 
size allow us to independently estimate the initial Lorentz factor of the jet. 

Subject headings: gamma rays: bursts 
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Introduction 



Gamma ray burst (GRB) afterglows are produced by the interaction of the ORB jet 
with the circumburst medium in which it is im mersed. An after glow spectrum is composed 



19981 ). and is produced when 



of a series of power laws separated by breaks (iSari et al. 
electrons in the circumburst medium spiral around the tangled and compressed magnetic 
field lines present at the shock between the jet and the circumburst medium. In this manner, 
the kinetic energy of the jet is gra dually converted into luminosity and particle energy 



( iMeszaros fc Rees 



1993 : 



Katz 



1994f ). The properties of the jet- including its initial kinetic 



energy, mass, and half-opening angle- only partially determine the temporal evolution 
of the afterglow. Because the afterglow luminosity is dependent upon the density of the 
circumburst medium, the nature of that medium (i.e., whether it is wind-like or a uniform 
density ISM) plays a crucial role in determining the tempora 



1998 



Huang et al. 



2000 



Chevalier fc Li 



2000 



Granot &: Sari 



evo 



2002 



ution as well ( 



Mesler et al. 



Sari et al. 



2012bf l. 



The flux density evolutio n of a GRB afte rglow is related to the density of the medium 



in which the GRB is present (jWaxman 



19971 ). The structure of the circumburst medium 



should be in some way imprinted upon a GRB afterglow, allowing the nature of the medium 
near the progenitor, and , henc e, the properties of the progenitor itself, to be identified 



( iMeszaros 



2002 



Woosley 



20111 ). Analytical models of relativistic jets in power law mediums 



(n(r 



oc r '^, with k = for a uniform density and k = 2 f or a wind) have yi e lded light 



curves that are in reasonable agreernent ^y ith observations (JBerger et al. 



2003 



Curran et al. 



2011 



Price et al 



20021 ) 



2000; 



Yost et al. 



Semi-analytic models developed by 



Frail et al. 



(120051) and 



Berger et al 



fl2003[ ) have 



been used to estimate parameters of the GRB 030329 radio afterglow, including the isotropic 
equivalent energy Ek, the electron spectral index p, the electron and magnetic field energy 
fractions (eg and e^, respectively), and the jet half-opening angle 9j, using broadband radio 
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observations over the course of the first year. 



van der Horst et al 



(120081) extend this to the 



first ~ 3 years. We use all previously published observations of GRB 030329 between 0.84 
and 250 GHz as well as archival VLA observations at 1.4, 4.9, and 8.5 GHz taken between 
1.7 and 5.5 years after the burst to perform burst calorimetry deep into the non-relativist ic 



regime. The broadband afterglow is modelled using a method developed in 



Mesler et al 



fl2012bh . 



2. GRB 030329 



At a redshift of 2; = 0.1685 ( Greiner et al. 



20031 ). the gamma-ray burst GRB 030329 



is one of the closest GRBs detected to date. Assuming a Acdm cosmology with Hq = 71 
km s~^ Mpc~^, Qm = 0.27, and ^a = 0.73, GRB 030329 is located at an angular distance 
of d-A = 587 Mpc, with 1.0 mas corresponding to 2.85 pc. The burst was first detected 
on March 29th, 2003 at 11:37 UTC by the High Energy Transient Explore 2 (HE TE- 2 ) 



satellite, and was subsequently localized in the optical by [Peterson fc Pried ( l2003l ). The 
GRB 030329 radio afterglow was (and still is) the most luminous afterglow to ever have 
been observed, achieving a maximum flux density of 55 mJy at 43 GHz one week after 
the burst. The relative proximity of the burst to the Earth, coupled with its extremely 



high luminosity, made it possible for the radio afterglow to 



Long Baseline Arra y (VLBA) at 5 GHz (jTaylor et al. 



3 e directly resolved by the Ve ry 



2004 



2005 



Pihlstrom et al. 



2007 



Mesler et al. 



2012al ). Additionally, the 5 GHz radio afterglow was detected by the VLA for 
more than 5.5 years. GRB 030329 has therefore provided a unique opportunity to study 
the evolution of both the luminosity and the physical size of a GRB afterglow 
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Archival Data 



We compare our model to observations of the GRB 030329 radio afterglow luminosity. 
In addition to previously published data, we include late-period 1.4, 4.9, and 8.5 GHz 
observations, taken from the NRAO data archive, that were made between 621 and 2032 
days after the burst (see Table [1]). The VLA data were reduced in the standard manner 
with the use of the AIPS software package. Either 3C138 or 3C286 was used for absolute 
flux calibration, and J1051+213 or J1021+219 were used to determine the complex antenna 
gains. 

To ensure the best possible fit of our semi-analytic model to the broadband afterglow, 
we also utilize previously-published measurements of the afterglow flux at 840 MHz, 1.4 
GHz, 2.3 GHz, 5.9 GHz, 8.5 GHz, 15 GHz, 23 GHz, 43 GHz, 100 GHz, and 250 GHz. 
Previ ously-published Wes t erbor k Sy nthesis Radio Telescope (^y SRT) data were taken 



from 



van der Horst et al. 



fl2005h and 



van der Horst et al. 



Very Large Array (V^ 



Pihlstrom et al 



A) d ata were taken from 



fl2007h . and 



Mesler et al. 



(VLBA) data point at 20 32 days comes from fJMesler et al. 



Berger et al 



(120081). P r eviously-pu b 



(120031) 



Frail et al 



ished 



(120051) 



(l2012ar). The 4.9 GHz V ery Long Baseline Array 



2012a). The 100 GHz and 250 



GHz data points are from 



Sheth et al. 



mom . 



4. Modeling the Jet Expansion 



In this section, we will discuss the semi-analytic method that we have de yeloped for 



modelli ng GRB afterglow emission. This model builds upon work appearing in 



Mesler et al. 



(l2012bh . It will be out 



since its publication in 



ined below with emp hasis on improvements that have been made 



Mesler et al 



( I2012bl ). 



gamma ray bursts are modelled as initially relativistic, double-sided jets that propagate 
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outward into an ambient, circumburst medium. In the following discussion, we will use 
unprimed quantities to refer to the reference frame that is comoving with the jet, we will 
use un-primed quantities with no subscript to refer to the reference frame in which the ISM 
is at rest, and we will use quantities with the subscript © to refer to the reference frame of 
an Earthbound observer. 

4.1. Jet Hydrodynamics 



An expression for the evolution of the jet Lorentz factor F can be found by invoking 



the requirement for the conservation of energy. The jet will sweep up materia 



circumburst medium as it propagates, forcing the jet to decelerate. 



Pe'er 



from the 



mm shows that 



dr 

dm 



7(r^-i)-(7-i)r/3^ 



M.. 



6771 



;i_e)^[2^r-(7-l)(l + r-2)] 



2^1' 



;i) 



where F is the Lorentz factor of the jet, Moj is the initial mass of the jet ejecta, 7n is the 
total mass that has been swept up by the jet, /3 = (1 — F~^) is th e normalized bulk 



velocity, 7 ^ (4F + 1)/(3F) is the adiabatic index (JHuang et al.ll2000l ). and 



syn 



is the radiative efficiency of the jet f lDai et al 



^ syn 



+ t'e 



(2) 



19991). The expansion time tL = t/(/3Fc 



and the synchrotron cooling time t' = Qnc/aT^e^T^pB (JDai fc Lu 



19981 ). where c is 



the speed of light, ar is the Thompson scattering cross-section, epsilo7ie is the fraction of 
the burst energy stored in the electrons, 7np is the proton mass, and B' is the magnitude 
of the comoving frame magnetic field. In the stellar wind medium case, the synchrotron 
cooling time is initially much larger than the expansion time, meaning that the jet starts 
off with a radiative efficiency e ~ e^. As the jet begin to propagate outward, the cooling 
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time eventually approaches and then surpasses the expansion time, driving the radiative 
efficiency to 0. The transition from partially radiative expansion (e ~ eg) to adiabatic 
expansion (e ~ 0) generally occurs after a few hours. The jet typically starts out adiabatic 
in the uniform density case and remains that way throughout its evolution. 

Our semi-analytic method is capable of producing light curves in arbitrary density 
profiles, but we limit ourselves for the sake of simplicity to a uniform density typical of an 
interstellar medium (ISM) and an n{r) oc r~^ stellar wind. Density profiles are modelled as 
a series of uniform grid points of width 10^"^ pc. The model is fed an initial kinetic energy 
Ek, ejecta mass Mej, and jet half-opening angle 6*0. The total mass swept up by the jet 
(m), the jet's Lorentz factor (F), the width of the leading edge of the jet perpendicular 
to the observer's line of sight (a), and the isotropic frame time (t) are then solved for 
simultaneously as a function of the radius (r) of the jet. 

The evolution of the jet is identical to the case where the GRB outfiow is isotropic until 
t = tjet, when the center of the jet comes into causal contact with its edge. For a jet with 
half-opening angle 6j = arctan a/r, this occurs when F ~ !/%• ^t. tjct, the jet experiences 
rapid lateral expansion, leading to an increase in the amount of circumburst material being 
swept up by the jet with time. The increase in swept-mass leads to a faster deceleration 
and a decrease in the jet luminosity. Eventually, the jet is decelerated to the point that it 
is no longer relativistic, and it transitions to spherical expansion. The transition time to 
non-relativistic expansion is denoted ^nr. 

4.2. Determining the Initial Lorentz Factor 



G RB 030329 is unique in t 



5 GHz flTavlor et al 



2004 



l at its afterglow 



2005 



Pihlstrom et al. 



las been directly resolved using VLBl at 

With multiple 



120071 : 



Mesler et al 



2012aD 



epochs of observation, the expansion history of the burst can be determined directly, 
allowing us to place constraints on the initial Lorentz factor of the burst. 

Initially, the mass being swept up by the GRB jet will be negligible as compared to the 
initial jet ejecta mass M^y The jet will coast at nearly constant speed in this regime. It is 
only when the two terms in the denominator of equation [T] are of the same order that the 
jet begins to decelerate appreciably. We will define the isotropic frame time t = tdcc as the 
time at which 

Mej = em + (1 - e) m (W) [27r - (7 - 1) (l + f-^)] . (3) 



The afterglow linear size evolves according to the comoving sound speed (JHuang et al. 



2OOOI) 



/7 (7-i)(r-i) 
i + 7(r-i) ■ ^ ^ 

In the regime where t < tdec, the Lorentz factor is nearly constant, meaning that the 
afterglow expands laterally at a nearly constant rate P±. Using the relationship between 
the isotropic frame time and the Earth frame time 

dt = T(T + Vr2 - 1) rfte, (5) 

we find that the relationship between the expansion rate P± of the afterglow and the initial 
Lorentz factor is 



^^^-^{ri^r.^) fK%'ir 
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for t < tdcc- The observed average lateral expansion rate (/3app) of the GRB 030329 radio 
afterglow is shown in Fig. [T^. At early times {t < tdcc), the lateral expansion rate is nearly 
constant, and, therefore (3± ~ (/^app)- After t = t^cc, however, the jet begins to decelerate 
and both f3± and (/3app) begin to decrease. The time tdec will therefore show up in the plot of 
(/^app) vs. t as a break where ^j"*" begins to decrease from its initial value of ~ 0. Looking 
at Fig. [1], it can be seen that (/^app) is consistent with a nearly constant value until after at 
least ~ 15 days, and begins to decrease at some point before 83 days. A transition to a 
slower rate of expansion appears to occur at some time between day 25 and 83. Therefore, 
we argue that 25 < tjcc ^ 80 days, with a fit to the data points taken at or before 83 days 
yielding P± = 3.9 ±1.9 (assuming P± is constant over this interval). Equation [6] can then 
be solved numerically to yield Fq = 3.4;*^]^'4. It should be noted, however, that it is possible 
to successfully fit a single power law to the data in Fig. [Dd. If the expansion rate were 
interpreted as decreasing monotonically over the entire time range, then our value of Fq 
becomes a lower limit on the initial Lorentz factor. It should also be noted that Equation [5] 
ignores the effects of time-dilation due to cosmological redshift, which is negligible for the 
nearby GRB 030329. 



The afterglow emissio n from GRB 



from a two-component jet (JBerger et al. 



330329 has previously been interpreted as coming 



2003 



Frail et al 



20051 ). One of the components was 



highly relativistic with a small initial half-opening angle of 6*0 — 5° and was responsible for 
the high energy emission (optical and higher frequencies), while the other component was 
only mildly relativistic with a larger jet half-opening angle 6q ^ 17°, and was responsible 
for emission at frequencies in the optical and below. Our value of Fq = 3.4l}4 is consistent 
with the Lorentz factor expected for the wider, moderately relativistic jet component. 

The initial Lorentz factor is related to the ejecta mass Mej and the kinetic energy Ek 
via 



<M 
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(a) 
i 

t 



10 10 10 

t (in days) 



t (in days) 



10 10 



(b) 

▼ 




\,^ T 



Fig. 1. — The observed size of the G RB 030329 radio afte rglow (a) and the observed trans- 



verse expansion rate (/3app) (b) from 



Mesler et al. 



( l2012al ). The sohd hne in b is composed 



of two fits. The horizontal component corresponds to the average value of the data points 
up to day 83, while the power law component corresponds to a fit of the data points from 
day 83 onward. 



Efc = (ro-l)Mejc2. 



(7) 



Determination of the initial Lorentz factor can therefore provide an important constraint 
on the ratio between the initial burst kinetic energy and the ejecta mass. 
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4.3. Fitting the Xr MHD Models 

Using the observations of the burst linear size (Fig. Wp), we produce a best-fit to our 
MHD models by varying the kinetic energy E^, the ejecta mass Mej, and the initial jet 
half-opening angle ^o- Best-fits to a variety of wind and uniform density profiles are shown 
in Fig. [21 We find that the afterglow size evolution can be successfully fitted to either a 
wind or a uniform density, but that lower values of Xr = X^ /^ ■, where A^ is the number 
of degrees of freedom present in the model, can be obtained for fits of uniform media to 
the afterglow size evolution than for stellar winds {xl uniform — ^ versus Xr wind — 8)- The 
medium surrounding GRB 030329 is therefore perhaps more characteristic of a uniform 
density ISM than a stellar wind. A good fit can be obtained for 2 < Fq < 10 over the entire 
range in n and A* that we searched (10~^ < n < 10^ cm~^ and 0.3 < A* < 6 cm~^). This is 
in agreement with the value of Fq = 3.4ti4 that we obtain using equation |6l 
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Fig. 2. — Fits of the GRB 030329 radio afterglow linear size to MHD models. Solid lines 
refer to models in which a uniform medium was assumed, while dashed lines indicate a stellar 
wind. 
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5. Modelling the Light Curve 

While the linear size evolution of the afterglow can give us insight into the initial 
Lorentz factor of the jet, it can only be used to determine the ratio between the initial 
kinetic energy and the ejecta mass if the nature of the circumburst medium is unknown. In 
order to determine the values of the kinetic energy and the ejecta mass, as well as the values 
of the other burst parameters, we must expand our model to include the propertie s of the 



Mesler et al. 



fl2012bl ) 



afterglow emission. We build upon the model detailed in section 5 of 

to model the GRB 030329 synchrotron emission by fitting to broadband radio afterglow 

observations between 840 MHz and 250 GHz taken between 1 and 2032 days after the burst. 

5.1. Spherical Emission and Beaming 

The leading edge of a GRB jet is not a fiat plane that propagates directly toward an 
Earthbound observer. Rather, it is s section of a spherical surface, meaning that material 
at different positions along the leading edge of the jet has different velocity components 
toward the observer. Relativistic beaming will cause radiation that is emitted by material 
with the largest velocity component in the direction of the observer to be beamed toward 
the observer more than radiation that is emitted by material with a slightly smaller velocity 
component. Additionally, jet material at different locations in the jet is not all at a uniform 
distance from the observer, leading to differences in light travel time throughout the jet. 
The overall effect of the relativistic beaming on the afterglow is an increase in the total 
observed flux at early times when the jet is still relativistic. The difference in arrival times 
of photons emitted from different regions of the jet makes the light curve broader and more 
smooth than it would be otherwise. 

To incorporate the effects of spherical emission and beaming into our model, equations 
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27-32 of iMesler et al.l (l2012bl ) must be modified to account for the fact that radiation 
emitted by different parts of the jet at the same time tcmit will not all reach the observer at 
the same time t^, nor will it all be beamed at the observer to the same degree. 

The arrival time t® at which radiation emitted by the jet reaches the observer will be 
dependent upon the angle 6 at which it was emitted with respect to a line connecting the 
center of the GRB progenitor to the observer. Material at angle 6 = will arrive at the 
observer first, while material emitted at 6* = 6j will arrive last. The equation for the equal 
arrival time surface is 



tfl 



1^1 — /3cosi 



dr = const. 



Numerical integration of equation [8] yields the radius at which emission at some angle 
6 was emitted in order to reach the observer at time t®. To account for the time and 
angular dependence of the beaming of radiation toward the observer, we integrate up the 
luminosities of the material at each location [r, 6, (p] which emit radiation that arrives at 
the observer at time t® to find the total afterglow flux density: 



F., 



47rL)2 






dcosO 



(9) 



where V is the Doppl er factor T> 



angle occupied by the jet ( iModerski et al 



1 — /3 cos 6) and Qj = 27r(l — cos 6j) is the solid 



20001 ). The quantity L'^^, is the co moving frame 



luminosity of the afterglow, which can be determined via equations 27-32 of 
mi2U \. 



Mesler et al. 
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6. Burst Calorimetry Prom the Broadband Afterglow 

The MHD and emission models detailed above were used to perform calorimetry on 
the broadband radio afterglow of GRB 030329. Seven parameters were fit simultaneously: 
the kinetic energy (Ek), the ejecta mass (Mej), the jet half-opening angle {Oj), the fraction 
of the total burst energy stored in magnetic fields (e^) and in electrons (eg), the electron 
power law index p, and the medium density (n) in the case of a uniform medium or the 
wind density scaling factor (A*) in the case of a wind. The factor A^, is set such that the 
medium density p(r) = 5 x lO^^A^r"^ g cm~^ for a stellar wind. 

The values that are obtained for the GRB 030329 burst parameters appear in Table |2j 
Our best fit models reproduce the overall evolution of the broadband afterglow with a few 
exceptions. At early times, scintillation effects produce fi uctuations in th e observed fiux. 



2005|) at 840 MHz and 



Large-scale refraction produces fiuctuations at late times (iFrail et al.l 

1.4 GHz. Taken together, these two effects, which are not accounted for in the emission 

model, are responsible for the large values of Xr that we obtain. 

Synthetic light curves produced using the best fit uniform medium and wind models are 
shown in Fig. [31 The best-fit uniform medium model predicts a medium density oi n = 7 
cm~^. The best-fit stellar wind model predicts a wind density scaling factor of A^, = 2.1. 

The best-fit uniform density and wind models are both capable of reproducing the 
broadband afterglow until about ^nr = 42 days, where the jet's transition to non-relativistic 
expansion produces a transition to a shallower decay in the fiux at all frequencies. The 
transition to non-relativistic expansion begins at roughly the same time in both the 
wind model and the uniform density model, but the length of time that it takes for the 
wind-model jet to become non-relativistic is approximately twice as long as the transition 
time in the uniform medium case. Hence, the decay in the wind-model fiux does not 
asymptote to the shallower power law as quickly as in the uniform medium case, and 
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the wind model cannot reproduce the real broadband afterglow after ^nr- We argue 



that a uniform density medium is suggested for the circumburst environment o 



030329, which is consistent with previous work (iBerger et al 



van der Horst et al. 



2005 



van der Horst et al. 



2008 



2003 



Mesler et al. 



Frail et al. 



2012a). 



GRB 



2005; 



Our fitted burst parameters are broac 



2003 



Frail et al 



2005 



van der Horst et al. 



y consistent with previous wor 



2005 



van der Horst et al 



i (IBerger et al. 



20081) (Table [3D. We 



find that the initial burst kinetic energy is Ek = 2.2 x 10^^ erg. Radiative losses and energy 
losses do to adiabatic expansion reduce this to Ek = 0.56 x 10^^ ergs by tjet = 12 days. 
The ejecta mass is found to be M^j = 2.4 x 10~^Mq, yielding Fq = 5.9. We also find that 
^0 = 27°, eB = 0.023, eg = 0.26, and p = 2.2. The value we obtain for the medium density of 
n = 7 cm~^ is somewhat larger than has been found by previous authors. The wide range in 
densities obtained for the GRB 030329 circumburst medium is probably due to differences 
in the individual models employed by each author. Uncertainties in the structure of the 
jet magnetic fields and the time dependence of e^ and e^ limit the accuracy of any gamma 
ray burst emission model. Given the imperfect nature of the current understanding of the 
jet physics, we estimate that the values we obtain for the burst parameters are accurate to 
within approximately a factor of 2. 

The uniform medium and wind models that produced the best fit to the broadband 
afterglow were themselves fit to the observed radio afterglow image size evolution (Fig. H]). 
Both models fit the evolution of the afterglow size, though the uniform medium fits slightly 



better (Xr, uniform = 6-1 vs xl 



wind 



7.3), again suggesting a uniform circumburst medium. 



From the initial kinetic energy and the ejecta mass, we find that Fq = 5.9 in the uniform 
medium case and Fq = 8.8 in the stellar wind case. The uniform medium model initial 
Lorentz factor agrees with our estimate of the initial Lorentz factor from section 2.2 of 
Fq = 3.4^J'4 to within ~ 1.6cr, while the wind model initial Lorentz factor agrees to within 
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3.4(7. 



Conclusions 



We have presented the first simultaneous fit of a gamma-ray burst afterglow model 
to a GRB's broadband light curve and its observed afterglow expansion rate. Up to 5.5 
years of afterglow observations at radio frequencies between 840 MHz and 250 GHz was 
employed to perform accurate calorimetry on the gamma-ray burst GRB 030329 deep into 
the non-relativistic phase. Values for the burst parameters were determined, and the nature 
of the circumburst medium was explored. 

By noting that a GRB jet will coast at a nearly constant velocity until a time tdec 
where it has begun to sweep up a significant amount of circumburst material, we derive 
a relationship between the average apparent afterglow lateral expansion rate, (/3app), and 



the initial Lorentz factor Fn. 



be Tfi 



3.4tii 



Berger et al. 



'he initial Lorentz factor for GRB 030329 is found to 



(120031 ) found that the GRB 030329 radio afterglow was 



best modelled as the product of a jet component that was initially mil dly-relat ivist ic . 



This interpretation 



van der Horst et al 



las been subsequently supporte d by other authors (JFrail et al 



2005 : 



van der Horst et al. 



2005 : 



20081 ). and is in agreement with our estimate 



of To. 

We find that a best fit of a model that assumes a uniform density circumburst medium 
to the broadband radio afterglow predicts values for the burst parameters that are similar to 
the values that are obtained when a stellar wind density profile is assumed. The goodness of 
fit Xr of the model that assumes a uniform density medium, however, is better by a factor 
of ~ 1.7 than the goodness of fit of the wind model, suggesting that a unifor m medium m ay 



be preferred over a stellar wind. This is in agreement with previous work f lBerger et al. 



-17- 



20031 : iFrail et al.ll2005l : iMesler et al.ll2012al ). but runs contrary to what is expected if GRBs 



have stellar progenitors. 

The best-fit wind and uniform density models were fitted to the observed average 
apparent afterglow lateral expansion rate and are both in good agreement with the data. 
The uniform density model agrees with our estimate of Fq = 3.4j;]^4 to within 1.6a, whereas 
the wind model agrees to within 3.4o", again perhaps suggesting that the circumburst 
medium is more characteristic of a uniform ISM than a stellar wind. 

The ability to determine the initial Lorentz factor of a GRB jet provides a powerful 
constraint on models of the burst's evolution. To determine the initial Lorentz factor, 
however, we must have enough observations of the afterglow image size to determine tdec- 
The GRB 030329 afterglow is currently the only gamma ray burst afterglow to have been 
directly resolved. In the future, more luminous, low-redshift bursts will need to be imaged 
with VLBl so that combined afterglow and lateral expansion evolution fitting can be applied 
beyond the case of GRB 030329. 

The National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 



-18- 



Date 


At 


Frequency 


Flux Density 


VLA Project Code 




(days) 


(GHz) 


(/iJy) 




2004 Dec 09 


621.2 


1.4 


650 ± 70 


AF414 


2004 Dec 09 


621.2 


8.5 


250 ± 30 


AF414 


2004 Dec 23 


635.1 


1.4 


590 ± 70 


AF414 


2004 Dec 23 


635.1 


4.9 


370 ± 60 


AF414 


2004 Dec 23 


635.1 


8.5 


300 ± 40 


AF414 


2005 Jan 23 


665.9 


1.4 


420 ± 80 


AK583 


2005 Jan 30 


672.9 


1.4 


480 ± 70 


AK583 


2005 Jan 31 


673.7 


1.4 


510 ±80 


AS796 


2005 Mar 31 


732.7 


8.5 


260 ± 40 


AK583 


2005 Apr 07 


739.7 


1.4 


650 ± 40 


AF414 


2005 Apr 07 


739.7 


8.5 


300 ± 30 


AF414 


2005 Jun 08 


801.5 


1.4 


380 ± 10 


AK583 


2005 Oct 21 


937.0 


8.5 


100 ± 30 


AK583 


2005 Oct 29 


944.5 


8.5 


120 ± 30 


AK583 


2005 Dec 14 


990.5 


8.5 


180 ± 20 


AK583 


2005 Dec 26 


1003.1 


8.5 


140 ± 30 


AK583 


2006 Mar 22 


1088.7 


1.4 


630 ± 30 


AS864 


2006 Mar 22 


1088.7 


8.5 


170 ± 20 


AS864 


2006 Apr 15 


1266.2 


8.5 


80 ±30 


AS864 


2006 Apr 29 


1126.5 


4.9 


080 ± 20 


AS933 


2006 Sep 15 


1266.2 


4.9 


170 ± 50 


AS864 



Table 1: Observations of the GRB 030329 radio afterglow at 1.4, 4.9, and 8.5 GHz taken 
from the NRAO data archive and not appearing in a previous publication. 
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Medium Type 


(10^1 erg) 


(Mo) 


n 


es 


ee 


p 


xl 


uniform (ISM) 


2.2 


2.4 X 10-^ 


27 


0.023 


0.26 


2.2 


309 


wind 


2.5 


1.8 X 10~'* 


21 


0.025 


0.38 


2.2 


526 



Table 2: Best fits to the GRB 030329 burst parameters assuming either a uniform density 
or a stellar wind circumburst medium. The data at all frequencies show significant scatter 
at early times due to scintillation effects. In addition, large-scale diffraction effects may 
be responsible for the scatter in the 840 MHz and 1.4 GHz data. These two processes are 
responsible for the large values of Xr- Note that the values of E^ listed here correspond to 
the kinetic energy of the jet at the time of the burst. 



Model 


(days) 


(days) 


Ek 
(10^1 erg) 


00 

n 


es 


ee 


P 


n 


Relativistic^ 


10 


N/A 


0.67 


26 


0.042 


0.19 


2.2 


3.0 


FulP 


14 


48 


0.90 


26 


0.074 


0.17 


2.2 


2.2 


Non-Relativistic^ 


N/A 


50 


0.78 


N/A 


0.13 


0.06 


2.2 


1.3 


Relativistic'^ 


10 


N/A 


0.24 


42 


0.43 


0.28 


2.2 


0.8 


Non- Relat ivist ic^ 


N/A 


80 


0.34 


N/A 


0.49 


0.25 


2.2 


0.8 


Full 


12 


42 


0.56 


27 


0.023 


0.26 


2.2 


7.0 



Table 3: Comparison of best-fit parameters for the various models that have b een produced 



for th e GRB 030329 radi o afterglow. These raodels we re first published in 



(120031)^ 



Frail et al. 



(120051 ) 



van der Horst et al. 



fl2005h 



van der Horst et al. 



Berger et al. 



fl2008f ll The 



bottommost model is the best-fit uniform density model presented in this work. Note that 
the kinetic energies listed here are valid at t = tjet for models valid in the relativistic regime 
and at t = ^nr for models that are valid solely in the non-relativistic regime. 
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t (in days) 




t (in days) 



Fig. 3. — Best fits to the GRB 030329 broadband radio afterglow assuming a uniform medium 
(a) and a stellar wind (b). 
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t (in days) 



Fig. 4. — Fits of the uniform and wind density models from Table [2]to the observed apparent 
lateral afterglow expansion rate (/3app). 
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